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Abstract 

The reaction pp/pp ^ tt +iet+X is an important background process for Higgs boson searches 
in the mass range below 200 GeV. Apart from that it is also an ideal laboratory for precision 
measurements in the top quark sector. Both applications require a solid theoretical prediction, 
which can be achieved only through a full next-to-leading order (NLO) calculation. In this work 
we describe the NLO computation of the subprocess gg ttg. 
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The reaction pp/pp /f+jet + A' is an important background process for Higgs boson searches in the mass range below 
200 GeV. Apart from that it is also an ideal laboratory for precision measurements in the top quark sector Both applications 
require a solid theoretical prediction, which can be achieved only through a full next-to-leading order (NLO) calculation. In 
this work we describe the NLO computation of the subprocess gg ttg. 



1. Introduction 

The main objective of the Large Hadron Collider 
(LHC) at CERN is the discovery of the Higgs bo- 
son and the measurement of its mass and couplings. 
To achieve this important goal a solid knowledge of 
the production mechanisms and the corresponding 
backgrounds is mandatory. In the Standard Model a 
light Higgs boson is currently favoured by the avail- 
able data. Using the recently updated top mass of 
m, = 178.0 ±4.3 GeV 1 1 1 the electroweak fits yield 
an upper bound of 251 GeV (at 95% C.L.) and a cen- 
tral value of niH = 117 GeV |2|. To achieve a high 
signal significance in the Higgs searches in general, 
different production and decay mechanisms are com- 
bined. In the range up to 200 GeV the so-called weak 
boson fusion (WBF) process with the subsequent de- 
cay of the Higgs into a W-boson pair plays a dominant 
role. The most important background for the WBF 
process comes from the ff +jet process |3|. A very 
precise knowledge of this process is thus mandatory 
for the discovery of the Higgs boson. It is obvious that 
for precise measurements of the couplings a precise 
background determination is equally important. For 
example, it has been shown in Ref . 1 4 1 that even if one 
assumes only a 10% uncertainty of the fF+jet cross 
section it is still the dominant theoretical uncertainty 
in the measurement of an = c?wbf y. B{H ^ WW). 
As also pointed out in Ref. 131 this accuracy might be 
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achievable only through a full next-to-leading order 
(NLO) calculation. In a recent analysis 15 161 , the pos- 
sibility to extract the background from extrapolation 
of experimental data has been studied. In this analy- 
sis it was found that a background determination with 
5-10% accuracy might be possible. This is a very 
promising result. On the other hand — given the sig- 
nificance of the precise background determination — 
we believe that a cross-check with a full NLO QCD 
prediction is important. At the very end — having 
a good understanding of both results — both meth- 
ods/results could and should be used as complemen- 
tary. 

In fact the ff + jet reaction is not only important as 
background for Higgs searches, it is also an impor- 
tant signal process on its own. It is well known that 
top quark physics allows a test of the Standard Model 
at high scale. In particular one can search for possi- 
ble extensions of the Standard Model at the scale of 
the top quark mass. As far as top quark production 
at hadron colliders is concerned, the state of the art 
is as follows. The differential cross section for top 
quark pair production is known to next-to-leading or- 
der accuracy in QCD [7' 8 '9' 10]. In addition the re- 
summation of logarithmically enhanced contributions 
has been studied in detail in Refs. I11I12I13I14I15I16I 
17'181. Recently also the spin correlations between 
top quark and antitop quark were calculated at NLO 
in QCD 1 191201 . Since single top quark production 
provides an excellent opportunity to test the charged- 
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current weak interaction of the top quark, it has also 
attracted a lot of interest in the past. In particular 
NLO corrections were studied in Refs. II21I22I23I24I . 
In Ref. [241 the NLO corrections for the fully differ- 
ential cross section are given, keeping also the spin 
information of the top quark. 

In that context the natural next step is the calcula- 
tion of the NLO corrections for ff +jet production. As 
far as top quark physics is concerned, interesting ob- 
servables to study are those that vanish if there is no 
additional jet. Such observables allow for a direct test 
of the dynamics in the top quark sector For example 
the asymmetry f25] 
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A^(f, forward) +A^(£, backward) 

is such an observable. Here A^(£, forward /backward) 
denotes the number of forward/backward-going lep- 
tons as a function of the rapidity y of the additional 
jet. The measurement of this asymmetry allows for a 
direct test of the production and decay mechanisms. 
Using similar observables one can search for exam- 
ple for anomalous top-gluon couplings. A more pre- 
cise understanding of the cross sections for pp/pp ~> 
tt+jets+X is also important for measurements of the 
top quark mass. 

2. Outline of the calculation 

In this section we briefly summarize the calculation 
of the NLO corrections for the subprocess gg — > tig. 
In view of the number of external legs and the top 
mass as additional parameter it is obvious that even 
partial results are in general quite lengthy. In the 
following we will restrict our attention only to those 
parts of the calculation where special care is needed 
to construct a numerically stable program. 

2.1. Virtual corrections 

The calculation of the virtual corrections proceeds 
via the following steps: 

1 . Generation of the Feynman diagrams using for 
example Feynarts |26| or QGRAF |27|. 

2. Reduction of the tensor integrals to scalar one- 
loop integrals. 

3. Reduction of the amplitudes to standard matrix 
elements. 



4. Numerical phase-space integration of the 
squared matrix elements, including appropriate 
phase-space cuts. 

Technically the most complicated part is the evalua- 
tion of the pentagon-diagrams. Two sample diagrams 
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Figure 1. Sample Feynman diagrams contributing to 
the virtual corrections. 



are shown in Fig.^ Let us first address the evaluation 
of the scalar 5-point integrals. To calculate these, we 
use two different methods. One calculation is based 
on the method given in Refs. 1281291 . The basic idea 
of this method is that finite 5-point integrals can be 
expressed in terms of 4-point integrals (see for exam- 
ple 13013 11321 '). To apply this observation also to soft- 
and mass-singular integrals they are rewritten accord- 
ing to Refs. 1 28 29 1 in the following way: 



-(mass.rf=4) _ p(mass,J=4) 
^sing. 



(2) 



Here denotes the original 5-point integral in d 
dimensions while £'(iTiass,c(=4) obtained from the 
original integral by dressing the massless propaga- 
tors with a small mass X. The subtraction term 
^sing^^ which has the same singular structure as 
the 5-point integral £'(mass,rf=4) -^^ jj^g YircAt A, 0, is 
obtained by studying the soft and collinear behaviour 
of £'(mass//=4) ^^jj expressed in terms of 3- 
point functions |33|. Rewriting now the finite inte- 
gral £'(mass,£/=4) jgj-jjjg 4-point integrals we thus 
succeeded in expressing the original 5-point integral 
in terms of 3- and 4-point functions. A more de- 
tailed discussion can be found in Ref. 1291. The sec- 
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ond method we used to calculate the five-point in- 
tegrals is based on the fact that, even for divergent 
integrals, it is possible to obtain a representation as 
linear combination of 4-point integrals (see for exam- 
ple Ref. 1 32 1). Expressing the 4-point function for 
t/ = 4 — 2e in terms of the finite 4-point function in 6 
dimensions plus a combination of 3-point integrals al- 
lows us also to shift all the divergences to the 3-point 
integrals. Defining the 5-point functions through 

E''ipo,pi,P2,P3,P4, m , '« 1 , m , mi , nn) 
\ [ ^ 1 

we obtain for example 

Eo{Q,PuPi ~ P3,P4- P2,-P2,m,,m,,Q,Q,m,)\^-^^ 

= P{tn)P{s45)Co{pi -p3,P4-P2,Pl,0,0,m,) 
+ P{t24)P{s35)Co{p4- P2,-P2,Pl -p3,0,OT,,0) 
- (fl3 - t24)^P{tn)P{t24)Pis35)P{s45) 

xCo{0,pi-p3,P4~P2,m,,0,0), (4) 

with P{x) = — mj) and Sjj — {pi + pj)^, tij — 
(Pi ~ PjY- The parton momenta are assigned accord- 
ing to g{pi)g{p2) t{p3)t{p4)g{p5). For the cases 
at hand it is possible to solve all the required box- 
integrals in 6 dimensions. We checked that the two 
methods yield the same results for the 5-point inte- 
grals E'' . 

Having solved the scalar integrals, the next step is 
the reduction of the 5-point tensor integrals to scalar 
one-loop integrals. In principle one could attack this 
problem using the standard Passarino-Veltman ap- 
proach |34|. This method leads to spurious singular- 
ities in individual terms at the phase-space boundary 
due to vanishing Gram determinants in the denomi- 
nator These spurious singularities create numerical 
instabilities when doing the phase-space integration. 
Note that the spurious singularities cancel if one com- 
bines the individual terms analytically before doing 
the numerical integration. One solution of this prob- 
lem is a time-consuming extrapolation technique, as 
was used for example in Ref. 1 29 1 . As an alternative to 
the extrapolation technique a different reduction pro- 
cedure 1351 was also used in Ref. 1291 . In this work 
we follow the method developed in Ref. f351. Es- 
sentially the same technique to reduce scalar 5-point 
integrals to scalar 4-point integrals is also applied to 



the tensor integrals. In this way the 5-point tensor 
integrals are directly reduced to 4-point tensor inte- 
grals. The explicit calculation shows that in this way 
the spurious singularities in individual terms, due to 
vanishing Gram determinants depending on 4 exter- 
nal momenta, are avoided |35 1. 

Let us just mention at the end that there are also 
other methods to solve the scalar 5-point integrals 
and perform the reduction of the tensor integrals. For 
example one could also use the methods developed 
in Refs. I36.37,.38J . (For Ref. |36J see also Walter 
Giele's talk in these proceedings.) 

2.2. Real corrections 

The calculation of the required matrix elements is 
straightforward. A sample diagram for the reaction 
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Figure 2. Sample Feynman diagram contributing to 
the real corrections. 



gg ttgg is shown in Fig. |2] We used two different 
methods to obtain the required colour-ordered helicity 
amplitudes: 

1 . A Feynman-diagram-based approach where we 
evaluate all the diagrams contributing to one 
specific colour-ordered subamplitude. 

2. Using the recurrence relations a la Berends and 
Giele |39|. 

We find complete agreement in the results of the two 
methods. Furthermore we also checked that our re- 
sults agree with the ones obtained using Madgraph 
1401 . To extract the singularities from collinear or 
soft partons we use the dipole subtraction method 
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r4r42'43V The idea of the subtraction method is to 
add and subtract a term which, on the one hand, can- 
cels pointwise the singularities of the matrix elements 
in the singular regions of the phase-space and is, on 
the other hand, easy enough to be integrated analyt- 
ically. Schematically the NLO contribution is then 
obtained from the following formula: 



C?NLO = 



m+1 



[cereal - OsubJ 









/ 


Ovirt. -f- 




J m 







finite 



dx / [afact.(^) +Csub(-^)] ■ 



(5) 



finite 



Here afact. {x) denotes the contribution from the fac- 
torization of initial-state singularities due to the pres- 
ence of coloured partons in the initial state. The con- 
tributions Osy|3,6sub obtained from a^ub by inte- 
grating out the 'unresolved' parton. The result is 
split into the two terms oJ^{,,asub to render the last 
two integrals individually finite. A remarkable fea- 
ture of the subtraction method is that the analytic in- 
tegration of the subtraction has to be done only once 
and that in the whole procedure no approximation is 
made. This is made possible by the universality of 
soft and collinear factorization in QCD. The explicit 
expressions for C^sub7 0^^j,, and Osub can be obtained 
from the colour-ordered subamplitudes using the for- 
mulae given in Ref . 1 43 1 . In particular Osub is obtained 
from a sum over individual dipole contributions. In 
the case at hand we have to include the contribution 
from 36 individual dipoles. We do not consider the 
splitting g ^ tt because the divergence is regulated 
by the quark masses. (For light quarks one could con- 
sider the corresponding dipoles to render the integra- 
tion numerically more stable.) We have checked that 
the combination of the 36 dipoles indeed reproduces 
all the singular limits arising from single unresolved 
configurations. 

3. Status and results 

The current status of the project is as follows. Most 
of the separate contributions are implemented in the 
form of computer programs allowing the numerical 
evaluation of the cross sections. In Fig.|3]we show as 
an example, the result at the parton level for the virtual 
corrections (defined as the second term in Eq. Q) for 
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Figure 3. Result for the virtual corrections for the 
subprocess gg ttg as defined by the second term in 
Eq. (ki_ > 20 GeV). 



different centre-of-mass energies. (Given that the sep- 
aration shown in Eq. (jSj involves some freedom, this 
individual contribution does not have a direct phys- 
ical interpretation unless the remaining contributions 
are added — we just show it for illustrative purposes.) 
As can be seen from Fig.|3lthe method we used for 
the treatment of the tensor integrals gives indeed nu- 
merically stable results. Furthermore we note that 
the inclusion of t/c'^j, together with the renormaliza- 
tion of the coupling and the quark mass renders the 
second term in Eq. ^ finite, as it must be. This is 
an important cross-check. As mentioned earlier we 
also checked that the integrand for the first contribu- 
tion in Eq. ^ is also finite for all single unresolved 
phase-space configurations. Given the complexity of 
the project we think it is very important to have in- 
dependent cross-checks for every individual contribu- 
tion. While most of the calculation is already cross- 
checked, we still work to finish also the remaining 
checks. Complete results will be presented elsewhere. 

4. Conclusions 

In this work we discuss the NLO calculation for 
the partonic reaction gg ttg. Up to remaining 
cross-checks, the calculation is almost finished. In 
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particular we have shown that the virtual corrections 
are stable using the reduction procedure discussed in 
Ref. |35|. 
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